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Effects of Charge Ordering on Spin- Wave of Quarter-Filled 
Spin-Density- Wave States 
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Spin-wave excitations of quarter-filled spin-density-wave state, which coexists with charge 
ordering, have been studied for one-dimensional extended Hubbard model with the nearest- 
neighbor repulsive interaction (V) and next-nearest-neighbor one (Via)- We calculate dispersion 
relations for the acoustic and optical spin-wave within the random phase approximation. Our 
numerical calculation shows that energy spectrum of the acoustic branch is well described by 
a simple sine function form. In the states coexisting with charge-density-wave, the spin-wave 
velocity decreases with increasing V or V2. Our numerical result, that velocity reduces to zero 
in the limit of large V or V2, is analyzed in terms of a spin 1/2 Heisenberg model with effective 
antiferromagnetic exchange interaction. 
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§1. Introduction 

Organic conductors, tetramethyltetraselenafulvalene 
(TMTSF) and tetramethyltetrathiafulvalene (TMTTF) 
salts,ElcP exhibit, due to an interplay of interaction and 
low dimension, exotic spin-density-wave (SDW) states, 
which coexist with the charge-density-wave (CDW), e.g., 
4fc F CDW in TMTTF salt and 2k F CDW in TMTSF 
salt (k F denotes a Fermi wave number). U n These states 
have been analyzed in terms of the extended Hubbard 
model with several repulsive interactions. Applying the 
mean-field theory, it is shown that the interactions with 
finite-range play a crucial role for the coexistence. The 
coexistence of SDW with 4k F CDW (2k F CDW) appears 
when the nearest-neighbor repulsive interaction V, (the 
next-nearest- neighbor repulsive interaction V 2 ) becomes 
large BcP 

The fluctuations around the SDW ground state have 
been examined extensively by calculating the collective 
modes within the random phase approximation (RPA). 
The modes for charge and spin fluctuations in the pres- 
ence of only on-site repulsive interactfon-l ID were evalu- 
ated for the incommensurate SDWQB&EiJ'H.^nd for the 
commensurate SDW state at quarter-fillmgjEj'EiJ where 
the latter may be relevant to above organic conductors. 
The commensurability energy induces an excitation gap 
in the charge excitation where the gap becomes zero at 
V = V C M> When V = a coexistent state of 2k F SDW 
and 4fcp CDW appears© and the harmonic potential for 
the charge fluctuation vanishes. The presence of V2 leads 
to an unusual ground state where 2k F SDW coexists jwi 
2k F CDW and also with Ak F SDW for V 2 > V 2c &\ 
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When V 2 = V 2c , one finds the disappearance of the har- 
monic potential with respect to the relative motion be- 
tween the density wave of the up spin and that of the 
down spin. For large V 2 , a new collective mode with 
an excitation gap appears where the mode describes the 
relative motion of density waves with opposite spin fol- 
lowed by amplitude mode_pf 2k F CDW. The excitation 
gap vanishes at V 2 = V^c-rJi 

While the spin-wave modes at quarter-fill iag: w as cal- 
culated for the conventional Hubbard model,E3Ef much 
of their problems is not known in the presence of V 
and V 2 . The excitation with long- wavelength has been 
calculated within the RPA for the case of V ^ and 
V 2 = 0,tZP and for the spin-wave velocity as a function 
of Vj_or V 2 v3 based on a functional integral formula- 
tionliS' The latter has shown that the spin- wave velocity 
decreases to zero (a finite value) with increasing V (V 2 ) 
when 2k F SDW coexists with Ak F CDW (2k F SDW coex- 
ists with 2k F CDW and 4/c F SDW) . In the path integral 
method, it remains an open question how tq_treat the 
coupling to long- wavelength spin fluctuations l3 

In the present paper, we examine the acoustic and 
optical spin-wave modes in the presence of V and V 2 
(leading to the charge ordering) for understanding the 
properties of the several kinds of transverse spin fluctu- 
ations of the SDW ordered states coexisting with CDW. 
We calculate not only the energy dispersion relation of 
them but also their collective operators within RPA. In 
§2, the mean-field ground state and the response function 
are explained for the calculation of the spin-wave modes. 
In §3, the resulting energy spectrum of them is given in 
conjunction with the characteristics of their collective co- 
ordinates including several spin fluctuation. In addition, 
the interaction dependence of the spin-wave velocities is 
shown as a function of V or V 2 . The effect of the charge 
ordering on the spin-wave is examined and is compared 
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with that obtained by the path integral methodic In 
§4, we discuss the effect of dimcrization and analyze the 
behaviors of the spin-wave velocity in the limit of large 
V or V 2 by using an effective Hamiltonian of a spin 1/2 
Heisenberg model. 

§2. Mean-Field Ground State and Response 
Function 

We study a one-dimensional extended Hubbard model 
at quarter-filling. The Hamiltonian is given by 

h = - 5> - {-iyt d ){clc J+1 . a + h.c.) 

+ H int , (2.1) 



N 

E 

3=1 



(Unj^riji + Vrijrij + i + l^rijUj+z) ,(2.2) 



where Cj CT denotes a creation operator of an electron at 
the j-th site with spin a = (f , |), and satisfies a periodic 
boundary condition C'j +Ar a — Cj a with the total number 
of lattice site AT, rij = + nji and nj a = Cj a Cj a . The 
quantity t and t d are the energy of the transfer integral 
and that of the dimerization. Quantities U, V and V 2 
correspond to coupling constants of repulsive interaction 
for the on-site, the nearest-neighbor site and the next- 
nearest-neighbor site, respectively. We take t and the 
lattice constant as unity. 

For a quarter-filled band, the mean-fields (MFs) of 
SDW and CDW are written as (m = 0, 1, 2 and 3, and 
sgn(a) = +(-) for a =\ (|)) 
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where Qo = 2fcp = tt/2 with fc F (= ?r/4) being the Fermi 
wave number. The z-axis is taken as the quantized axis. 
The expression < O ^mf denotes an average of O over 
the MF HamiltonianJiEP given by 



Hmf = 
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Fig. 1. Phase diagram in the V-V 2 plane for U = 4 and t d = oil 
Three kinds of states correspond to a pure state of 2/cp SDW(I), 
a state of 2fc F SDW and 4fc F CDW (II) and a state of 2k F SDW, 
2fcp CDW and 4fcp SDW (III). These states are shown schemat- 
ically in the region where the arrow (circle) represents the spin 
(charge) density at each lattice site. The dashed curve denotes 
a first-order transition between regions II and III. 



where e k = — 2£cosfc. The expressions of each MF are 
given such that S = 0, D = 1/2, S Qo = S , 3 * Qo = Sie ie , 
D Qo = D* Qo ee £ lC i(e ^ /2) , S 2Qo = S* 2Qo ee S 2 and 
D 2Qo = D 2 Qo = D 2- Quantities Si(> 0), D x {> 0), 
S 2 and D 2 denote amplitudes for 2fc F SDW, 2£; F CDW, 
4fc F SDW and 4fc F CDW, respectively. Their values de- 
pend on the corresponding-ground states. The quantity 
9 denotes a phase of SDW.H 

The ground state has been obtained previous!} 
where the phase diagram on the plane of V 2 and V is 
shown in Fig. 1. When V 2 = and t d — 0, a pure 
2fc F SDW state (6> = tt/4, = tt/2) (I) is found for 
V < V c (~ 0.34) and a coexistent state of 2k F SDW^md 
4fc F CDW {6 = 0, (j> = tt/2) (II) is found for V > 
where cj) = tan" 1 (S'i/£)i). When V = and t d = 0, the 
pure 2/c F SDW state (I) changes to a coexistent state of 
2fc F SDW, 2fc F CDW and 4i F SDW (6 = tt/4, < tt/2) 
(III) for V 2 > V 2c (~ 1.32)0 Note that the finite-range 
interactions V and V2 contribute to the appearance of 
the quantities D\ and D 2 . In the following, we use the 
states I, II and III for those in the regions I, II and III, 
respectively. 

The spin- wave is examined for the above three ground 
states. The spectrum of the spin-wave modes is calcu- 
lated from the pole of the response function given byfcf 



-1 

27V 



/3 

dT(T T $(q,T)&(q))> 



(2.5) 



where w„(= 27m//3) and (3~ 1 {= kB.T) denote the 
Matsubara frequency and temperature, respectively. 
The components of operator &(q) consist of the 
transverse spin fluctuations given by S am (q) = 
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J2-ir<k<n^k ej ^k+q+ m Q {a = x,y) where ip\ = 
(Cfct)C^.) and a a (a — x,y) are Pauli matrices. For 
the state I (0 = 7r/4) and the state II (0 <-B < 7r/4 due 
to id), the operator &(q) can be taken aaiiJ 



s; 2 (9)) 



(2.6) 



We note that, for small g, two components S^g) and 
^3 (?) are related directly to the transverse spin fluctua- 
tion of the SDW with wave numbers Qo and — Qq respec- 
tively and that other components of 5^ (g) and Sy 2 (q) 
denote the fluctuations correlated with those of the wave 
numbers and 2Qp. The coupling to long- wavelength 



spin fluctuation Sy (q) gives a 
spin- wave velocity multiplied by a factor [1 — fW(O)] 1 / 2 
in the weak coupling limitEf where iV(0) is the density of 
state per spin at Fermi level. For the coexistent state of 
2fc F SDW, 2k F CDW and 4/c F SDW (the state III) with 
the spin modulation with Qo and 2Qo, one needs to ex- 
tend the operator $ T (g) into that with eight components 
written as 



renormalization of the 



^(q) = (St 1 (q),SUq),St (q): 
^1(9), ^3(9), 5^ (g), 



St 2 (q), 



j y 3yrU' u yow;i ^y2\<u) ■ (2-7) 

In addition to S xl (q) and S^ 3 (g), the component S^ 2 (g) 
for small g is also related directly to the transverse spin 
fluctuation of the SDW with 2Qo while other compo- 
nents comes from coupling to these three components. 
Operators of eqs. (2.6) and (2.7) act to create the trans- 
verse spin fluctuations by raising or lowering the spins 
within the same sites (i.e., rotation of single electron due 
to on-site interaction). However, the finite-range inter- 
actions induce the transverse spin fluctuations with the 
spin current (i.e., spin rotation followed by neighboring 
electrons) , which ia^related to the bond spin-density- wave 
(BSDW) state.lEfEy In the present paper, we neglect the 
effect of spin current, which is left for the study in the 
separated paper together with the problem of BSDW. 
Based on such a consideration, eq. (2.5) within RPA is 
calculated as 



^RPA 

R (g,w) = 



n (q,u) 



1 + U n (q,ui) 



(2.8) 



where the polarization function, n (q,uj), is evaluated in 
terms of the MF Hamiltonian, cq. (2.4). 



n (q,u) = 



-1 
2N 



dr(T T $(g,r)$t( ? )) 



MF 



x e iw " r |. ... . (2.9) 

The pole of eq. (2.8) determines the excitation spectrum 
of the collective mode, w 7 (g), i.e., 



det(l + C/n (g,w 7 (g))) =0 , 



(2.10) 



where 7 = Tl or T2 and ujti {^T2) denotes the acoustic 
(optical) mode. The eigenvector, $ 7 (g), which describes 
the collective mode, is obtained by 



l + [/n(g,w 7 (g))U> 7 (g) = 



(2.11) 



Actually, in terms of ^(g) and $ 7 (g), the operator for 
the spin- wave with w 7 (g) is written as 



77 7 (g) = $ 7t (g)$(g) 



(2.12) 



an explicit form of which will be given in_the next section. 
For u) ~ w 7 (g), eq. (2.8) is rewritten ast3 



A-> (g,w 7 (g)) 



2 w 2 - W 2( g ) ' 
2 adj (l + UU (g,w 7 (g)) 
(det(l + J7n (q,w)) 



S( _d 



(2.13) 



(2.14) 



where the matrix A 1 (g,w 7 (g)) is proportional to the 
residue and adj {l + U Yl (g,w 7 (g))^ is the adjoint ma- 
trix. The spectral weight of the collective mode is exam- 
ined by calculating following diagonal element, 



Al(q) 



(?,"(?)) 



/U 



(to = 1,3,0,2) . 

(2.15) 



When td = and V = V2 =0, one finds that, in the limit 
of small U, Af 1 ^) = 4P(0) = 2tto f with v F = y/2 and 
AP(0)=AP(0)=A2(0)=A2(0) = 0. 

We examine the V (or V2) dependence of the velocity 
of acoustic spin-wave given by 

uri(q) 



v = lim 

«— 



q 



(2.16) 



which is compared with the previous work 







§3. Spin- Wave Modes 

The characteristic properties of the spin-wave modes 
are examined for the various kinds of SDW ground states, 
which are given by the states I, II and III. The numerical 
calculation is performed by taking U — 4 and t — 1 where 
id = in this section. 

3.1 Effect of V 

First, we calculate the spin- wave spectrum w 7 (g) as 
a function of V for V2 = 0. In Fig. 2, the spectrum of 
both the acoustic and optical spin-wave modes obtained 
from eq. (2.10) are shown for V — 0, 0.5, 1.5 and 2.0 
with V% — and id = where the state II is obtained for 
V > V c (— 0.34). Our numerical calculation within the 
visible scale shows that spectrum of the acoustic mode 
takes a simple sine function form with respect to the wave 
number g, i.e., WTi(g) = ^>ti{Qo/2) sin(2g). The optical 
mode, WT2(g), is well separated from the continuum for 
the parameters in Fig. 2. The spectrum of the state II is 
similar to that of the state I. 

The spectral weights for the acoustic mode, A^(q), 
are shown in Fig. 3 where Af 1 (q) and A^ 1 (q) (AQ 1 (q) 
and A^ 1 (q)) are weights for spin fluctuations of x (y) 
direction with wave numbers Qo + q and —Qo + g (g 
and 2Q + g) , respectively. For the spectral weights with 
q ~ or Qo, the dominant weights are given by the flue- 
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Fig. 2. Excitation spectrum of the spin-wave modes, u>Tl(g) 
(acoustic mode) and u)T2{l) (optical mode) for U = 4, V2 = 
and td = 0. 




Fig. 3. Normalized spectral weight of the acoustic mode, 
A^ l 1 (q)/(2TTV F ), for V = 0, 1.5 and 2.0 in Fig. 2. 



tuations with ~ ±Qo (note that S y o(Qa) = S y i(0) and 
S y 2(Qa) = ^3(0)), as shown for a model with only UM> 
Figure 3 indicates that the deviation of the spin orien- 
tation arising from the spin-wave changes continuously 
from the x direction to the y direction when q increases 
from to Qq. It is found that all the weights are sup- 
pressed by V. 

Here we examine the behavior in the limit of small q. 
From the numerical result of eq. (2.11), we found that 
the eigenvector $ T1 (q) for the acoustic mode is written as 
$ T1 (<? -> 0) oc (e ie ,e- ie ,0,0) where 9 = tt/4 for V < V c 
and 8 — for V > V c . Then, the spin- wave operator of 
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Fig. 4. Schematics of spin-wave for q — > are shown, (a) acoustic 
mode in the state I, (b) acoustic mode in the state II, (c) optical 
mode in the state I and (d) acoustic mode in the state III. The 
lower part of each figure represents a projection of the spin on 
the x-y plane. 



eq. (2.12) in the limit of q — > may be expressed as 

VT1 (0)<xY,S*cos(Q r ] +0) , (3.1) 



ip\(T a ^j/2 and ip\ = (Cj^C^). Equation 
(3.1) denotes that the gapless excitation at q = is de- 



where S" 
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Fig. 5. V dependence of the gap, 1^7^2(0), for U = 4, V2 = and 
td = where u>T2(q) disappears at the location shown by the 
arrow. The dotted curve denotes the bottom of the continuum. 
In the inset, the corresponding spectral weight, (Q) / (2irv-p) 
is shown and the weight is zero if not shown explicitly. 



Fig. 6. V dependence of the spin-wave velocity v, for U = 4, V2 = 
and td = where the dot-dashed curve-denotes the previous 
result obtained by a path integral method]!^ In the upper inset, 
v° as a function of U while, (0) / (2nvp) corresponding to the 
main figure is shown in the lower inset. 



scribed by rotating the spin quantization axis toward 
x direction uniformly. The spatial dependence is illus- 
trated in Figs. 4(a) and 4(b) for V < V c and for V > V c , 
respectively, which are written simply as (/* /* ,/ f) f° r 
V < V c and {■/■/) for V > V c . Similarly, w T i(Qo), 
corresponds to the excitation with a uniform rotation 
toward y direction. 

Another spectrum 10x2(1) of the optical mode has a 
gap. The corresponding eigenvector $ T2 (q) is given by 
& T2 (q — > 0) oc (e w , -e~ w , 0, c) where 9 = ir/4 and c > 
for V < V c , and 9 = and c = for V > V c . Thus, it is 
found that the spin-wave operator for the optical mode 
in the limit of q — > is given by 

77 T2 (0) oc [ 2 <SJ sin(Q r i + 9) + icS] cos(2Q r J )] . 

3 

(3.2) 

For V < V c , the spatial dependence is shown schemati- 
cally as where a plane of (\/*) is not parallel 
to that of (\ /) due to c 7^ 0. The explicit spin config- 
uration is illustrated in Fig. 4(c). For V > V c , it is given 
by (<— J— »J,) where the symbol <— denotes a component 
of Sj. In Fig. 5, the gap of the optical mode, u>T2(0), 
is shown as a function of V. The thin dotted curve de- 
notes the bottom of the continuum. With increasing V, 
wt2(0) moves into the continuum at the location shown 
by the arrow, The inset shows the corresponding spectral 
weight. The spectral weight, (m=l,3) reduces and 
becomes zero at V 1 corresponding to the arrow, i.e., 
above which the optical mode disappears. For q — > 0, 
the optical mode in the state I contains the fluctuation 
of S V 2(q) in addition to the fluctuations of S x ±(q) and 
S X 3{q) while S y 2(q) is absent for the mode in the state II. 

Now, we examine the spin- wave velocity v, which is 
evaluated by using eq. (2.16). Figure 6 represents the 



V dependence of v where the dot-dashed curve denotes 
the result obtained by the path integral methodJMJ The 
quantity v , which stands for the velocity in the case of 
U = 4 and V = V 2 = t d = 0, is given by v° = 0.39 (1.29) 
for the solid (dot-dashed) curve. In the upper inset, the 
U dependence of v° is shown, where the solid curve and 
dot-dashed curve correspond to those in the main figure. 
There is a small jump at V — V c due to the discon- 
tinuity of Si and Z?2, which occurs only for id = 0E 1 
With increasing V(> V c ), v (solid curve) takes a slight 
hump and decreases to zero more rapidly than that of the 
dot-dashed curve. The coupling to long-wavelength spin 
fluctuations was neglected in the path integral method. 
Therefore, it is found that the velocity is strongly renor- 
malized by the coupling to such a spin fluctuation. The 
hump is not essential since it disappears for t<j ^ (see 
Fig. 10). The behavior of v for large V can be understood 
as that of spin 1/2 Heisenberg model with the effective 
exchange interaction J ~ t A /(UV 2 ) : which is discussed 
in §4. The lower inset shows the corresponding weights 
of A^(0) (m = 1,3), which decreases with increasing V 
in a way similar to v. 

Here we comment on the fact that the spin-wave ve- 
locity as a function of V changes rapidly for V > V c 
while the velocity remains the same as that of V = for 

V < V c . The latter behavior of the velocity originates 
in the present scheme of treating the response function. 
For the case of V < V c , in which the pure 2kp SDW 
state (I) is realized and the MF is independent of the 
nearest-neighbor interaction V, the quantities n (q,oS) 
m eq. (2.9) and R m (q ,u>) in eq. (2.8) do not depend on 
V. If we consider, however, the transverse spin fluctu- 
ations with the spin current which come from the term 
of the nearest-neighbor interaction V, it is expected that 
the corresponding response function R (q, lo) explicitly 
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Fig. 7. Excitation spectrum of the spin- wave mode, wti (q), with 
some choices of V% for U = 4, V = and t<j = 0. 



Fig. 8. q dependence of the elements of eigenvector | 9 ri (q) |, cor- 
responding to the case of V2 = 2.0 in Fig. 7. 



includes the interaction V though the extended polariza- 
tion function n (<Z,w) is still independent of V. Thus, 
the V dependence of the spin-wave velocity may appear 
in the state I. 

3.2 Effect of V 2 

Next, we calculate the spin- wave mode as a function of 
V 2 in the state III by using eqs. (2.7) and (2.10). Figure 
7 shows the excitation spectrum of the acoustic mode 
wti (q) with some choices of V 2 for U — 4, V — and = 
0, where the coexistent state of 2kp SDW, 2/cp CDW 
and Ak F SDW is found for V 2 > V 2c (^ 1.32). The 
optical mode uiT 2 {q) moves into the continuum leading 
to a disappearance of the mode for V 2 ( ^ 1.35) being just 
above V 2c . With increasing V 2 (> V 2c ), the spectrum 
u>Ti{q) decreases. We note that the dispersions for V 2 = 
1.5 and 2.0 is also well described by a relation, Wti(?) = 
w T i(Qo/2)sin(2g). 

The corresponding eigenvector <3? T1 (g) for V 2 = 2.0 is 
shown in Fig. 8. For the calculation of the eigenvec- 
tor in cq. (2.11), we choose seven components except for 
Syi(q) in $(9) of eq. (2.7) since there are two degener- 
ate modes, i.e., the mode with dominant contribution 
from S xm (q) and the mode with that from S yrn (q). The 
element |$?i(<?)| corresponding to S y i(q) is negligibly 
small, and then the spectrum calculated from seven com- 
ponents is actually nearly the same as that of Fig. 7. It is 
worthy to note that the mode with small q in the region 
III of Fig. 1 is composed of not only S xm (q) (to = 1,3) 
but also S x2 (q) (Fig. 8). Such a feature can be under- 
stood by considering the spatial configuration, in which 
the spin modulation is of type (tjJ,J.) (Fig. 1) and the 
variation of spin moment comes from the presence of 
4/cf SDW. In this case, the spin-wave operator express- 
ing the gapless excitation at q — > is expected to take a 



form of 

Vti(0) = Y. S J t 2Bl cos (Qorj +6) + B 2 cos(2Q r j )] , 
3 

(3.3) 

which means a uniform rotation toward x direction with 
keeping the relative angle between the each spin, like 
(/ y /)■ Here, the coefficients B\ and B 2 satisfies 
the relation B 2 jB\ = S 2 / S\ according to the spin mod- 
ulation in the ground state. In eq. (3.3), the term of 
B 2 comes from the fluctuation of S x2 {q). Actually, our 
numerical calculation shows that 



Tl 



oc 



(e ie , e" ie , 0, c', 0, 0, 0) , (3.4) 



q^0 



where 9 = tt/4 and c' = |$i 2 1 (0)|/|$^ 1 1 (0)| = 5 2 /5i (~ 
1.2 for V 2 — 2.0). The spatial dependence of this spin- 
wave is explicitly illustrated in Fig. 4(d). 

The spin-wave velocity v corresponding to that in 
Fig. 7 is shown as a function of V 2 in Fig. 9 (solid curve). 
The dot-dashed curve denotes the result obtained pre- 
viously by the path integral method. With increasing 
V 2 (> V 2c ), the solid curve decreases to zero rapidly while 
the dot-dashed curve tends to a finite value (v — > t) for 
large V 2 S3 Thus the present RPA calculation indicates 
the importance of the couplings to long-wavelength spin 
fluctuations, which has been neglected in the previous 
path integral method. The behavior in the limit of large 
V 2 will be discussed in the next section. 

§4. Discussions 

We have examined spin-wave modes in the presence 
of charge ordering, which exhibits a clear effect of sup- 
pressing the spin-wave velocity even without dimeriza- 
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Fig. 9. V2 dependence of v, for U = 4, V = and id = 0. The Fig. 10. V dependence of v, for U = 4, V$ = and the fixed 



dot-dashed curva-denotes the previous result obtained by a path 
integral method.1 1 8 P 



id = 0, 0.1, 0.3 and 0.5. The dotted curves represent v = 2Jj = 
2aj[(i 2 — id 2 ) 2 /(UV 2 ), where an ~ 1.6 is estimated from the 
inset. 



tion (t<j = 0). Since the dimerization enhances the ef- 
fect of the repulsive interaction, one may expect fur- 
ther suppression of the spin- wave velocity by adding the 
dimerization. In this section, we discuss such an effect of 
dimerization (id 7^ 0) together with the limiting behavior 
of large V or V2 ■ 

First, we plot the velocity v as a function of V in 
Fig. 10, with the fixed i d = 0, 0.1, 0.3 and 0.5 (U = 4 and 
V2 = 0). The dimerization id has an effect of reducing 
the velocity, which is qualitatively, consistent with the 
result of the path integral meth.od.Ef The V dependence 
of v in the region II of Fig. 1 can be fitted by a dotted 
curve, which is proportional to (i 2 - t 2 d )/(UV 2 ). We 
analyze such a case by noting that the ground state with 
the spin and charge modulation given by (0, f , 0, J., 0) and 
(0,1,0,1) already shows up for even if id 7^ 0. 

In this state II, the antiferromagnetic exchange coupling 
between the up and down spins is induced by the process 
of the fourth order of the transfers t ± id leading the 
effective Hamiltonian as 



m 



N/2 

Jn ^2 S 2 j ■ S 2 j- 



l 1 f— 





O.'l 




/^"o 


:/ 


0.3 " 


' ■ ■ ■ ■ 




0.5 1 

2 2 2 

h)i(uvi) . 



(4.1a) 



Fig. 11. V2 dependence of v, for U = 4, V = and the fi xed id 
0, 0.1, 0.3 and 0.5. The dotted curves represent v = 2 yj Jm : 
1.6 X 2(t 2 — td 2 ) 2 /(£/V 2 2 ), which is estimated from the inset. 



2 \2 



Jn 



Mt 2 - g ) 
uv 2 



(4.1b) 



with Sf 



where a n = 12 and Sj = (Sf, Sj, 5| 
tp^(j a if!j/2 (ipj = (Ct|, Cj|)). This expression of Jn with 
an = 6 has been obtained by Mori and Yonemitsu.0- ) 
From eq. (4.1a), the spin- wave dispersion is given by 
u>(q) = Jn sin(2o) and then the spin-wave velocity is 
given by v = 2Jj. By comparing v = 2Jj with the 
velocity calculated from eq. (2.16) in the inset of Fig. 10, 
we found an ~ 1.6. Such a value of an differs apprecia- 
bly from that of the effective Hamiltonian although one 
obtains a good, coincidence between them for the case 
at half-filling .E3 It is not yet clear at present why it is 



complicated to estimate quantitatively by the effective 
Hamiltonian at quarter-filling. 

Next, we examine V2 dependence of v with some 
choices of id, which is shown in Fig. 11. With increas- 
ing id, v is suppressed monotonically in the region III 
(Fig. 1). For arbitrary value of id in the region III, v can 
be well fitted by a dotted curve which is proportional to 
(i 2 — id) 2 /(C^^2 2 )- Such a behavior can be also analyzed 
in terms of the following effective Hamiltonian. In this 
case, the spin and charge modulation in the ground state 
is given by (0, |, I, 0, 0, T, I, 0) and (0, 1, 1, 0, 0, 1, 1, 0), re- 
spectively. The exchange coupling between the spins of 
(|, I) mainly comes from the second process of the trans- 
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for, while that between the spins of (j, 0, 0, |) comes from 
the sixth order process. Thus, the effective Hamiltonian 
for large V 2 can be expressed as the Heisenberg model 
with the alternating interactions, i.e., 



N/4 

-ffJn* ~ E {JmSij-2 ■ Sij-i + J^Sij-i ■ S4.7+2) 

3=1 



am{t - t d ) s 
U 



, _ a' m (t + t d ) 2 (t 2 -t 2 ) 2 



p 

■ 1 if — 



(4.2a) 
(4.2b) 

(4.2c) 



where an = 4 and a' w = 1. In this case, one can 
find the expression of the spin-wave dispersion (Ap- 
pendix) as , uj(q) = v /J m J^sin(2g), where y/Jm J£ = 
^a w a' m {t 2 - tl) 2 /(UVf). By comparing v = 2y/JmJ^ 
with the numerical result in the inset of Fig. 11, we find 
x/amaJn — 1-6, which is compatible with that of the ef- 
fective Hamiltonian, i.e., 2. 

In summary, we have calculated the spin-wave excita- 
tions for a one-dimensional model with a quarter-filled 
band to investigate the effects of the nearest (V") and 
next-nearest- neighbor (V 2 ) interactions on them. From 
the dispersion relations for spin-wave modes in the coex- 
istent state of SDW and CDW, we found that the spec- 
trum of acoustic mode and then the spin-wave velocity 
are reduced by CDW (induced by V or V 2 ) and that the 
spin-wave in the coexistent state for large V or V 2 could 
be described by spin 1/2 antiferromagnetic Heisenberg 
models with the effective exchange coupling. 
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Appendix: Spin- Wave of the Heisenberg Anti- 
ferromagnets 

By taking four lattice sites as a unit cell and intro- 
ducing the sublattice A and B, the Hamiltonian of eqs. 
(4.1a) and (4.2a) can be rewritten as 



off 



N/4 

E 

1=1 



(JSb,i ■ Sa,i + J Sb,i ■ Sa,i+i) , (AT) 



where J = J' = Jn for eq. (4.1a), and J = Jn and 
J' = J m for cq. (4.2a). By using the Holstein-Primakoff 
transformation, 





= 1/2-0101 : 






(A-2) 


$A,l 


= S Al + \S y Al 


= (1- 


a\ai) 1/2 ai , 


(A-3) 


$A,l 


_ QX -QV 

— ^A,l — ia A,l 


= 4(i 


-a\ ai )^ 2 , 


(A-4) 




= -1/2 + ^6/ 






(A-5) 




= Sb,i + iS B> i 


= bj(l 


-^) 1/2 , 


(A-6) 



SB^S^-iSl^il-blh) 1 ^ , (A-7) 

cq. (AT) up to the second order with respect to a q and 
b q is given by 

+ (J + J'e- i4 ")6t a t + ( j + J' c ^)b q a q ] , (A-8) 

where the constant term is subtracted and the boson 
operators a q and b q are given by 

a\ = (N/4)-^ 2 J2e iqRl a\ , (A-9) 
1 

b\ = (N/A)- 1 ' 2 c~ iqRl b] ■ (A- 10) 

1 

By applying the Bogoliubov transformation to eq.(A-8), 
the spin- wave dispersion, u>(q), is obtained from 

u(q)-(J + J')/2 -(J + J'e i4 «)/2 \ 
( J + J'c- i4 «)/2 uj(q) + ( J + J')/2 ) - U ' 

(ATI) 



det 



Equation cq. (ATI) leads to the spin- wave dispersion for 
the Hamiltonian (AT) as 

Lu(q) = VjJ 7 sin(2q) , (AT2) 

where ui(q) = Jnsin(2o) for the Hamiltonian (4.1a) and 
w(o) = V JmJ^sm(2q) for the Hamiltonian (4.2a). 
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